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ABSTRACT 

The first short Chandra and XMM-Newton observations of the young and energetic pulsar J 1357-6429 pro- 
vided indications of a tail-like pulsar-wind nebula associated with this object, as well as pulsations of its X-ray 
flux with a pulsed fraction pf > 50% and a thermal component dominating at energies E < 2 ke V. The elon- 
gated nebula is very compact in size and might be interpreted as evidence for a pulsar jet. The thermal radiation 
is most plausibly emitted from the entire neutron star surface ofa 10 km radius and a 1.0±0.1 MK temperature, 
covered with a hydrogen atmosphere. At higher energies the pulsar's emission is of a nonthermal origin, with 
a power-law spectrum of a photon index T = 1.1 ± 0.2. This makes the properties of PSR J 1357-6429 very 
similar to those of the young pulsars Jl 1 19-6127 and Vela with detected thermal radiation. 
Subject headings: pulsars: individual (PSR J 1357-6229) - stars: neutron — stars: X-rays 
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1. INTRODUCTION 

Once discovered in radio, young rotation-powered pulsars 
with characteristic ages of r c = P/2P ~ 10 kyr (P is pulsar spin 
period) often become targets for X-ray observations because 
these objects usually possess large values of the rotational- 
energy loss, E > 10 36 ergs s" 1 , and hence are expected to gen- 
erate observable nonthermal (magnetospheric) emission and 
power detectable pulsar-wind nebulae (PWNe; see Kargalt- 
sev, Pavlov, & Garmire 2007 for properties of a number of 
PWNe). Measuring parameters of the magnetospheric radia- 
tion is particularly important for inferring emission processes 
operating in pulsar magnetosphere. Studying a PWN is essen- 
tial for determining the energy of the pulsar wind, probing the 
ambient medium, and understanding the shock acceleration 
mechanism(s). 

On the other hand, neutron stars (NSs), born hot in super- 
nova explosions, cool down via neutrino emission from the 
entire NS body and heat transport through the stellar enve- 
lope to the surface and subsequent thermal radiation. The NS 
cooling rate is determined by various processes of the neutrino 
emission, as well as the (unknown) properties of the NS inner 
matter (Yakovlev & Pethick 2004). Till recently, young pul- 
sars were thought to be such powerful nonthermal X-ray emit- 
ters that their thermal radiation would be completely buried 
under the nonthermal component. However, new observa- 
tions revealed the pulsars Vela (r c = 1 1 .3 kyr) and Jl 1 19-6127 
(t c =1 .6 kyr) as examples of young NSs with the bulk of the X- 
ray flux of a thermal origin (Pavlov et al. 2001; Gonzalez et 
al. 2005). The surface emission is of a special interest because 
confronting observational data with model thermal radiation, 
presumably formed in NS atmospheres, can allow one to in- 
fer the surface temperatures and magnetic fields, constrain the 
NS mass (M) and radius (R), and understand the properties of 
the superdense matter in the NS interior (see Zavlin 2007a 
[hereafter Z07] for a review). 

At a distance of d ~ 2.5 kpc derived from the pulsar dis- 
persion measure and the Galactic electron density model 
of Cordes & Lazio (2002), PSR J 1357-6429 (P = 166 ms, 
t c = 7.3 kyr, £ = 3.1 x 10 36 ergs s" 1 , surface magnetic field 
B~8x 10 12 G), discovered by Camilo et al. (2004; here- 
after C04) during the Parkes multibeam survey of the Galac- 
tic plane, is one of the nearest objects in the group of about 



two dozen currently known pulsars 1 with r c < 15 kyr. The 
actual pulsar's age could be r = 2r c /[n- 1] < 15 kyr, for a 
typical value of the magneto-dipole braking index n ~ 2-3. 
Radio data on this pulsar (C04) revealed a strong glitch of 
its spin period, AP/P ~ -2.4 x 10~ 6 , similar to those expe- 
rienced by some other young pulsars (e. g., Vela), support- 
ing the hypothesis that PSR J 1357-6429 is a young NS. It is 
located relatively close to the supernova remnant (SNR) can- 
didate G309.8-2.6, but a possible connection between these 
two objects has yet to be established (see C04 for details), op- 
posite to the case of the majority of young pulsars associated 
with SNRs. 

As with many other young pulsars, PSR J 1357-6429 has 
been been a target of both Chandra and XMM-Newton. Sec- 
tions 2 and 3 describe these observations and analysis of the 
collected X-ray data. The results are discussed in Sec. 4. 

2. OBSERVATIONS AND DATA REDUCTION 

XMM-Newton observed PSR J 1357-6429 on 2005 August 
5 for 11.6 and 14.5 ks of effective exposures with the Euro- 
pean Photon Imaging Cameras, EPIC-pn and MOS (respec- 
tively), operated in Full-Frame Window mode. The high- 
est achieved time resolution of 73.4 ms is not sufficient for 
timing of a signal of a 166-ms period. The data were re- 
processed with the XMM-Newton Science Analysis Software 
(SAS v. 7.0.0). Due to strong particle flares a background 
level during about 85% of the total observational span was 
higher by a factor of 3-20 than the "quiescent" rate. 

The pulsar was observed with the Chandra High Resolution 
Camera (HRC-S) operated in Timing mode on 2005 Novem- 
ber 18 and 19, for 16.1 and 17.1 ks of effective exposure 
and with a 16-^s time resolution. The Chandra Interactive 
Analysis of Observations (CI AO) software (v. 3.4; CALDB 
ver.3.3.0.1) was used to generate the level 2 event files from 
the two observations. Those were combined with the CI AO 
script , merge_all' , that resulted in the data product of 
a 33.2 ks exposure. The pulse-invariant (PI) channels with 
PI<25 or PI>150 were filtered out to minimize background 
contamination. 

3. OBSERVATIONAL RESULTS 
' jhttp : / /www . atnf . csiro . au/research/pulsar| 
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FIG. 1 .— Chandra HRC-S 16" X 16" image smoothed with a 0"4 FWHM 
Gaussian showing PSR J1357-6429 surrounded by an extended structure 
elongated in the east-north direction (indicated by the white-dashed box of 
a 1" X 1"5 size). White contours correspond to intensity values of 0.17, 
0.55, and 1.74 counts arcsec -2 ks~'. 



3.1. Spatial analysis 

Figure 1 presents a smoothed HRC-S image of the 16" x 
16" region around the radio position of PSR J 1357-6429 de- 
rived by C04. The image reveals a point-like source sur- 
rounded by a weaker diffuse emission. The source is centered 
at the position R.A.= 13 h 57 m 2?54 and Dec.= -64°29"30f'0. 
The difference of 0."74 between the X-ray and radio posi- 
tions is very close to the la error in the Chandra HRC-S po- 
sitional astrometry 2 . The diffuse emission feature, extending 
toward northeast to a distance of about 2" from the center of 
the point-like source, is seen in the both HRC-S datasets and 
is not an image artifact caused by the instrument anomalies 
and/or background 3 . The presence of the extended emission 
is also apparent in Fig. 2 showing two 1-D distributions of 
counts in rectangular apertures shifted along the main symme- 
try axis of the box plotted in Fig. 1 : one histogram presents the 
distribution of detected counts, the other indicates the HRC-S 
point-spread function computed from a 2-D simulation per- 
formed with the CI AO ' mkpsf ' tool. The proximity of the 
X-ray and radio positions and the morphology of the extended 
feature strongly suggest that the detected X-ray radiation is 
emitted by PSR J 1357-6429 and its possible PWN 4 . 

The 1"05 x 1 ''58 box shown in Fig. 1 contains 18 counts, 
of them only 14% belongs to background. Hence, the elon- 
gated "tail" of the diffuse emission (within the box) emits 
at a rate of 0.48 ±0.14 counts ks" 1 . Modeling the tail's 
emission with a power-law (PL) spectrum of a photon in- 
dex r = 1.5 (typical for X-ray PWNe — e. g., Kargaltsev et 
al. 2007) absorbed with the hydrogen column density «h,2i = 
«h/(10 21 cm" 2 ) = 5.0 (see Sec. 3.3) yields the tails's (unab- 
solved) flux F t ~ 3.3 x 10" 14 ergs cm" 2 s _1 in 0.5-10 keV, or 
the luminosity L t ~ 2.5 x 10 31 iif 5 ergs s _1 (^2.5 = <i/[2.5kpc]). 

No large-scale signatures of a diffuse emission which could 




2 |http : //asc . harvard. edu/cal/ASPECT/celmon| 

3 |http : //cxc . harvard. edu/proposer/POG/html/HRC . htnil| 

4 No other point-like source as a possible counterpart of this feature is 
found in available observational catalogues. 



FIG. 2. — Distribution of counts detected with the Chandra HRC-S in- 
strument (solid histogram) in rectangular apertures with a 1"84 height and 
varying widths (0."13 or 0"40) shifted along the main symmetry axis of the 
box shown in Fig. 1 . The dashed histogram shows the HRC-S point-spread 
function computed for the same apertures. 



be associated with the SNR candidate G309. 8-2.6 (Sec. 1) 
were found in these data. An upper limit on the averaged 
surface brightness of the possible SNR emission, estimated 
from the 2f5-radius circle centered at the position R.A.= 
13 h 56 m 57?0 and Dec.= -64°28"00, is 5 counts ks -1 arcmin" 2 . 

3.2. Timing analysis 

For the timing analysis, 137 counts were extracted from 
the l."05-radius circle centered at the pulsars X-ray position, 
with 132 counts estimated to belong to the pulsar. The photon 
times of arrival were corrected to the solar system barycenter 
using the CI AO ' axBary' tool and the pulsars radio posi- 
tion and convolved with the accurate timing solution obtained 
from radio observations of the pulsar conducted at the Parkes 
telescope 5 in the interval from 2005 April 27 to 2005 Decem- 
ber 10 (53,487-53,714 MJD): / = 6.019298216(5) Hz and, 
/ = -1.29968(6) x 10" 11 Hz s" 1 at the epoch 53,694.0 MJD 
(given with the lcr uncertainties; F. Camilo 2007, private com- 
munication). This radio solution leaves a timing residual of 
only 3.6 ms (2.2% of the spin period). At this ephemeris, 
the Z 2 -test (m is number of harmonics) results in the value 
of Z 2 = 12.0 (Z 2 = 14.6), or a signal detection at a 99.752% 
confidence level. 

Figure 3 presents the pulse profile of PSR J 1357-6429 ex- 
tracted at the radio spin parameters, with the estimated pulsed 
fraction, pj = 63 ± 15% (defined as the fraction of counts 
above the minimum in the light curve). The x 2 -test yields the 
value x 2 = 4.07 for v = 9 d.o.f, meaning a 95.996% probabil- 
ity that the pulse profile differs from a constant line. Because 
of the poor statistics, the shape of the pulse profile is not well 
determined. The only rather certain fact is that the pulsed 
fraction is large, implying special properties of the detected 
X-ray emission (see Sec. 4). 

3.3. Spectral analysis 

To optimize the signal-to-noise ratio, the EPIC spectra were 
extracted from 20"-radius circles centered at the pulsar posi- 
tion (these include the emission of the faint PWN suggested 
by the HRC-S data). The extraction regions contain about 

5 The Parkes Observatory is part of the Australia Telescope, which is 
funded by the Commonwealth of Australia for operation as a National Facil- 
ity managed by CSIRO. 



PSR J 1357-6429 in X-rays 
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FIG. 3.— Pulse profile of PSR J 1357-6429 extracted from the Chandra 
HRC-S data. The dashed line indicates the mean level. 
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FIG. 4.— Spectra of PSR J1357-6429 detected with the XMM-Newton 
EPIC instruments and fitted with a two-component, NS atmosphere (NSA) 
plus power law (PL), model. 



75% of the source flux. The instrumental responses were gen- 
erated with the SAS x rmfgen' and "arfgen' tools. Back- 
ground was evaluated from similar regions in vicinity of the 
pulsar position. The estimated source count rates (corrected 
for the finite extraction aperture) are 56 ± 5 and 17 ± 2 counts 
ks _1 in 0.3-10 keV, for EPIC-pn and MOS, respectively. 

A single PL model fits the spectra rather well, with min xt = 
1.17 for v = 59 d.o.f, yielding T = 2.1 ±0.4 and « H ,2i = 2.3 ± 
0.8 (lcr errors are given here and below). The derived non- 
thermal (isotropic) luminosity in 0.5-10 keV is L non ± ~ 2.0 x 
10 32 d 2 5 ergs s _1 = 6.5 x 10~ 5 Edl 5 (of this, about 10% belongs 
to the PWN). On the other hand, it is reasonable to assume 
that PSR J 1357-6429 emits observable thermal radiation. 
Adding a thermal component improves the fit, minx 2 =0.76 
for ^=57 d.o.f, indicating that it is required by the data at a 



99.9998% (4.8<r) confidence level. The thermal component 
can be equally well described with both a blackbody (BB) 
spectrum and a NS magnetized (B = 1 x 10 13 G) hydrogen at- 
mosphere model 6 . However, as in many other cases, the esti- 
mates on the temperature and size of the emitting area derived 
from these two models are very different(see Zavlin 2007a 
for details): = 1.1 ± 0.2 MK and R™ = (2.5 ± 0.5)^2.5 
km (redshifted values) in the BB fit, and T eff = 1 .0 ± 0. 1 MK 
(effective temperature) for R = 10 km at o?2.5 = 1 in the NS 
atmosphere fit. The corresponding bolometric luminosities 
are Lj£ ~ 3.7 x 10 32 <ff 5 ergs s" 1 and L atm ~ 7.1 x 10 32 ergs 
s" 1 (L~ = [l-2GM/c 2 /?]L atm ~ 0.6L atm forM= 1.4M and 
R = 10 km). The parameters of the nonthermal component 
are similar, V = 1.3 ±0.2 and 1.1 ±0.2, in the BB and at- 
mosphere fits (respectively), with L non ± ~ 1.4 x 10 32 c/2 5 er g s 
s"'= 4.5 x \Qr 5 Ed\ 5 . The inferred hydrogen column density, 
«H,2i = 4.9 ±2.0, is somewhat larger than in the single-PL in- 
terpretation but still consistent with the standard estimate sug- 
gested by the pulsar dispersion measure DM = 127.2 cm" 3 pc, 
nH,DM/(10 21 cm" 2 ) = DM/(10 20 cm- 2 ) ~ 4.0 (assuming a 10% 
ionization degree of the interstellar hydrogen), and is reason- 
ably lower than the galactic neutral hydrogen density esti- 
mated in the direction to PSR J1357-6429, n H i/(10 21 cm" 2 ) ~ 
10.0. It is also worth mentioning that the thermal componet 
contributes 72% of the X-ray flux observed at E < 2 keV, 
regardless of whichever model is used to describe the com- 
ponent. The best two-component fit involving the NS atmo- 
sphere model is shown in Fig. 4. 

4. DISCUSSION 

The first X-ray observations of PSR J 1357-6429 have re- 
vealed interesting properties of this object 7 , including the in- 
dication of the compact tail-like PWN associated with the pul- 
sar. The estimated efficiency with which PSR J 1357-6429 
powers the PWN, L t /E k, 0.8 x 10" 5 , can be considered as 
"medium" among those measured for tail-like PWNe associ- 
ated with other pulsars. For example, it is lower by about an 
order of magnitude than the efficiencies of the PWNe of the 
young pulsar B 1757-24 (Kaspi et al. 2001) and the millisec- 
ond pulsar B 1957+20 (Stappers et al. 2003; Zavlin 2007b) 
but larger by a factor of a few than those inferred for the 
Geminga's tail (Pavlov, Sanwal, & Zavlin 2006) and the 
Vela's southest jet (Pavlov et al. 2003). On the other hand, the 
derived luminosity of the PWN associated with PSR J 1357- 
6429 is about by a factor of 5-8 smaller than the nonthermal 
luminosity emitted by the pulsar itself, that does not match 
the correlation, Lpwn ~ 5L nont h, inferred by Kargaltsev et 
al. (2007) for a sample of pulsars. It is premature to make firm 
conclusions about the properties of the possible PWN. How- 
ever, despite the fact that the proper motion of PSR J 1357- 
6429 has yet not been measured, the shape and estimated lu- 
minosity of this elongated feature may suggest a hypothesis 
that it is a pulsar jet along the pulsar spin axis aligned with its 
velocity, similar to the interpretation proposed for the tail-like 
features of the Crab, Vela, and Geminga pulsars (Weisskopf et 
al. 2000; Pavlov et al. 2003; Pavlov et al. 2006). On the other 
hand, if the elongation direction and that of the proper motion 
turn out to be perpendicular to each other, it would indicate 



The 'nsa' code in the XSPEC package. 
7 Esposito et al. (2007) analyzed the same data, obtained similar results 
in the PL and BB-plus-PL spectral fits and put an upper limit pf < 30% for 
a possible sinusoidal signal. They searched for a diffusion emission of a 10" 
size and failed to notice the small extended feature. 
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FIG. 5. — NS cooling models (Yakovlev &Pethick 2004) with the '2p' pro- 
ton superfluidity (solid curves) and without the superfluidity (dashed curves), 
for different NS masses. The crosses indicate the pulsars: Jl 119-6127 (1), 
J1357-6429 (2), Vela (3), B 1706-44 (4), J0538+2817 (5), and B2334+61 (6). 
The redshifted surface temperature is 7;°° = [1 -2GM/c 2 R] 1 / 2 T eff (= 0.777/ eff 
for M/Mq = 1.4 and R = 10 km), with T e ff being the effective temperature of 
the NS surface obtained in the fits with the NS atmosphere models. 



that the X-ray emission of the PWN emerges from a torus as- 
sociated with the equatorial pulsar wind. Future measurement 
of the pulsar proper motion is important for validating this hy- 
potheses. 

The spectral data on PSR J 1357-6429 show that the bulk 
of the pulsar X-ray flux is probably of a thermal origin. It 
makes this object the second youngest NS in the group of en- 
ergetic pulsars with a thermal component dominating at E < 2 
keV. In addition to PSRs Jl 1 19-6127, J 1357-6429 and Vela, 
the group also includes PSRs B 1706-44 with r c = 17.5 kyr, 
J0538+2817 with the derived true age r ~ r c /20.5 = 30 kyr 
(Kramer et al. 2003), and B2334+61 with t c = 40.9 kyr (see 
Z07 for details). Similar to the situation with other members 
of this group, the origin of thermal component of PSR J 1357- 
6429 cannot be unambiguously determined from the spectral 
data alone. The parameters of the thermal component ob- 
tained in the BB fit may be interpreted as radiation from a 
small hot area (polar caps) on the NS surface, although the es- 
timated radius, ~ 2.5 km, is much greater than the canon- 
ical polar cap estimate, /? pc = (2nR 3 / cP) l l 2 ~ 0.3 km. Oppo- 
site to this, the fit with the magnetized NS atmosphere models 
indicate that the thermal emission originates from the entire 
NS surface. One could discriminate between these two inter- 
pretations invoking pulsations of the X-ray emission. In the 

8 This assumption should be taken with caution — see the case of 
PSR J0538+2817. 



case of PSR J 1357-6429, as the Chandra HRC-S instrument 
is most efficient at E < 2 keV, the pulsed fraction of p/ > 50% 
suggests that the thermal emission is intrinsically anisotropic, 
as predicted by the NS atmosphere models — otherwise the 
effect of strong gravitational bending of photon trajectories 
near the NS surface would strongly suppress the pulsations. 
This obviously contradicts the simplistic (isotropic) BB inter- 
pretation. Such a pulsed fraction also indicates that the NS has 
a strong nonuniformity of the surface temperature and mag- 
netic field. Hence, the temperature r e ff inferred from the NS 
atmosphere fit assuming the uniform surface should be con- 
sidered as an approximate estimate on the "mean" surface 
temperature. The same conclusions stand for PSR Jill 9— 
6127 with a pulsed fraction p j > 60% (Gonzalez et al. 2005), 
as well as PSR 0538+2817 whose pulse profile of the ther- 
mal flux emission shows a complicate dependence on photon 
energy (Zavlin & Pavlov 2004). 

Assuming that the characteristic ages of the five pulsars of 
this group are close to the true ones 8 , it is illustrative to com- 
pare NS cooling models with the surface temperatures esti- 
mated for these objects from the NS atmosphere fits. Figure 5 
presents cooling models with and without proton ('2p' model) 
superfluidity in the NS core (Yakovlev & Pethick 2004), as 
well as the measured temperatures. The superfluidity reduces 
the neutrino emission by suppressing the Urea processes and, 
hence, decelerates the NS cooling. This effect depends on NS 
mass, being stronger for higher masses. The comparison sug- 
gests that the interiors of these six pulsars are superfluid, and 
the NS masses may be in the range M = [1.5-1.6]M© (note 
that this range would be different for another model superflu- 
idity). 

However, despite the important X-ray information on 
PSR J 1357-6429 available, the pulsar's properties remain yet 
rather uncertain because of the scanty number of photons de- 
tected in these short exposures. Therefore, deeper observa- 
tions are required to provide data of much better quality to 
firmly establish the actual shape, size, and spectrum of the 
possible PWN, and accurately infer the X-ray spectral and 
temporal properties of the pulsar. 
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